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We report a study of the decays B° — > J/ipK ( -*^°ir + n~ , which involve the creation of a uu or dd 
quark pair in addition to a b — > c(cs) decay. The data sample consists of 110 pb _1 of pp collisions 
at yfs — 1.8 TeV collected by the CDF detector at the Fermilab Tevatron collider during 1992-1995. 
We measure the branching ratios to be BR(B° -> J /ipK*°-K + -k') = (8.0 ± 2.2 ± 1.5) x 10~ 4 and 
BR(_B° -> J/ipK°ir + TT~) = (1.1 ± 0.4 ± 0.2) x 10~ 3 . Contributions to these decays are seen from 
^{2S)K { * )0 , J/i/)K°p , J/iPK'+tt-, and J/ipKi(1270). 

PACS numbers: 13.25.Hw, 14.40.Nd 



The measured inclusive branching ratio for B — > J/ipX 
of (1.16 ± 0.10)% is considerably larger than the sum of 
the individual branching ratios of the known exclusively 
reconstructed decays [Q . One possible source of the miss- 
ing decay modes is a class of decays in which a quark pair 
is created in addition to a b — ► c(cs) decay. The CLEO 
collaboration recently reported observing one such mode, 
B — > J/tpepK 0, which involves an ss quark pair. This 
analysis studies similar B° decays that involve uu or dd 
quark pairs, an example of which is shown in Fig. m. 




FIG. 1: Example of a 6 — > c(cs) decay with a dd quark pair 
and spectator. 

These modes are potentially useful for CP violation 
measurements. For example, J/ipK s p° is accessible from 
both B° and B° which allows CP violation due to inter- 
ference between decays with and without mixing. An 
angular analysis of the decay products could separate 
the CP-even and -odd contributions. If the factorization 
hypothesis H holds, a time dependent analysis of the 
angular correlations could measure cos(2/3) in a manner 
similar to B° -> J/yj(K*° -> K^ir ) @. 

This analysis was performed using pp collisions 
recorded with the CDF detector, which is described in 
detail elsewhere |5| . For this analysis the important com- 
ponents are the Silicon Vertex Detector (SVX) , the Cen- 
tral Tracking Chamber (CTC), and the central muon 
systems. The SVX provides a track impact parameter 
resolution of ~ (13 -I- 40/px) /xm, where pr (in GeV/c) 
is the component of the momentum transverse to the 
pp collision axis (the z axis) Q. The CTC is a drift 
chamber whose charged particle momentum resolution is 
Spr/Pr ~ 0.001/(GeV/c). Two muon systems separated 



by 60 cm of steel cover the region |ry| < 0.6 for muons 
with pt > 1.4 GeV/c. Each of these central muon sys- 
tems consists of four layers of planar drift chambers. The 
inner system is separated from the interaction point by an 
average of five interaction lengths of material. An exten- 
sion to the central muon systems covers 0.6 < \rj\ < 1.0. 

This analysis uses a three-level dimuon trigger. The 
first level selects events with two separate sets of at 
least three linked hits in the muon chambers that cover 
1 77 1 < 1.0. The second level requires drift chamber tracks 
with px > 2.0 GeV/c which extrapolate to the linked 
hits in the muon chambers. The third level accepts 
J/tp — > fi + \i~ candidates with a reconstructed invari- 
ant mass between 2.8 and 3.4 GeV/c 2 . In addition to 
this trigger path, approximately 10% of the events which 
pass the level 3 dimuon trigger come from single muon 
triggers at levels 1 and 2 with muon px thresholds of 
either 7.5 GeV/c or 12 GeV/c. 

The offline analysis reconstructs J/ip — > p + pT candi- 
dates but does not require them to be the same candidate 
which passed the trigger. The offline reconstruction only 
uses muons which intersect both of the muon systems 
that cover \r}\ < 0.6. 

A ratio of branching ratios is measured between a sig- 
nal mode B° -> J/if:K^ Q 7T + 7r and a well established 
reference mode B° -> J/yjK^ 0: 

BRjjjg _ £re/ N s j g 
BR, re y ^sig N re f 

Many systematic uncertainties cancel in this ratio. The 
ratio of efficiencies R e = tref/^sig is determined with a 
Monte Carlo simulation. The number of signal and ref- 
erence events (N S i g and N re t) are measured in the data 
while applying similar selection criteria to both signal 
and reference decay modes. The only selection criteria 
which differ are those placed upon the two extra pions 
of the signal events for which there are no equivalents in 
the reference decays. 

From the dimuon trigger events, B° decay candidates 
are selected which satisfy the basic topology of the de- 
cays of interest. The reference modes are reconstructed 
by combining a J/ip — » candidate with either a 

Kg — » 7r + 7r~ or K*° — > K + ir~ candidate. K s candidates 
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are required to have an invariant mass between 485 and 
510 MeV/c 2 and to point back to the J/tp decay vertex. 
Additionally, the K s candidates' decay vertices are re- 
quired to have a positive displacement in the xy plane 
from the J/tp decay vertex with at least 5<r significance. 
K*° candidates are required to have an invariant mass 
between 820 and 970 MeV/c 2 and to originate from the 
J/tp candidate decay vertex. To reduce backgrounds, all 
final state particles are required to havepr > 0.5 GeV/c, 
to be within AR = yj ' Aif + A<j) 2 < 1.0 of each other, 
and to originate within 5 cm of each other in z. The sig- 
nal modes are reconstructed identically to the reference 
modes, except for the addition of two pions which origi- 
nate from the J/tp candidate decay vertex. The invariant 
mass of these two extra pions, m(7T + 7r~), is required to 
be greater than 0.55 GeV/c 2 . 

The final vertex fit constrains the particles to originate 
from a common vertex, except the Kg daughters whose 
combined momentum must point back to that vertex. 
For this fit, the invariant masses of the K s and J/tp can- 
didates are constrained to their world average measured 
values ji| . The \ 2 of the fit is required to have a confi- 
dence level above 0.1%. 

If there are multiple B decay candidates in the same 
event, all candidates are kept. Multiple candidates have 
the largest effect in the B° — > J/tpK*°ir + ir~ sample 
where misassignments of the K and tt from the K*° result 
in a broad Gaussian shaped background. Multiple candi- 
dates do not significantly affect the B° — > J / tpK^ir+n^ 
sample. 

To reduce background levels, additional selection cri- 
teria are placed on the transverse momentum of the 
neutral kaon, px(K^*'°), the proper decay time of the 
B° candidate, ct(B), and a B° isolation variable, I = 
Pt(B)/(pt(B)+pt(x)). The quantity px (x) is the scalar 
sum of the transverse momenta of all non-£?° candidate 
tracks within AR < 1.0 of the B° candidate momentum 
direction. For each signal mode these selection criteria 
are optimized to maximize S 2 /(S + Bkg) of the signal 
sample where S is the expected signal size and Bkg is 
the expected background size. This optimization uses 
the sidebands of the invariant mass distribution of the 
signal data and the invariant mass distribution of the 
reference data. 

The optimized selection criteria for B — > 
J/tpK 



- are p T (K*°) > 2.4 GeV/c, 
ct{B) > 170 (im, and I > 0.60. Figure |(a) shows 
the resulting invariant mass peak. The data are fit using 
the sum of a narrow and a broad Gaussian and a linear 
background. The width of the narrow Gaussian is fixed 
to 11.8 MeV/c 2 , based upon the expected width from 
the Monte Carlo simulation scaled up by the ratio of 
widths between the data and Monte Carlo simulation 
for the reference mode B° — > J/tpK*°. The fit results 
in 36.3 ± 9.9 signal events. The B° -> J/tpK*° refer- 
ence mode with similar selection criteria has 257 ± 18 



signal events. The ratio of reconstruction efficiencies 
between B° -> J/tPK*° and B° -> J/ipK* <K + ic- is 
R £ = 3.75. These numbers lead to a ratio of branch- 
ing ratios of BR(B° -> J/tPK* Q ir+ir-)/BR(B° -> 
J/ipK*°) = 0.53 ± 0.15, where the error is statisti- 
cal only. Using BR(B° -> J/tPK*°) = 1.5 x 10~ 3 
jlj, this corresponds to a branching ratio of 
BR(J/^*°7r+7r-) = (8.0 ± 2.2) x 10~ 4 . 




5 5.1 5.2 5.3 5.4 5.5 5.6 

JA|/K , °7t + jrMass [GeV/c 2 ] 



5 5.1 5.2 5.3 5.4 5.5 5.6 
JA|/ Ks Jt + 71" Mass [GeV/c 2 ] 



FIG. 2: Invariant masses of B° candidates for 
J/ipK*°w+-K- (a) and B° -> J/tpK^ir+ir- (b). 



B° 



The final state J/^jK*°tt + tt could come from the in- 
termediate states B° tP(2S)K*° or B° -> J/tPK*°p°. 
Within ±2cr of the B° invariant mass shown in Fig. ||(a), 
there are 9 tp(2S) — > J/ipir + ir~ candidates on an ex- 
pected background of 3 within ±10 MeV/c 2 of the ^(25) 
mass of 3.686 GeV/c 2 [0. The remaining signal can- 
didate events have higher J/tpTr + Tr~ invariant masses. 
There is no identifiable resonant structure in the tt + tt~ 
invariant mass distribution to indicate a large p° contri- 
bution, nor is there any identifiable resonant structure to 
the K *°7r ± or K* tt + it~ invariant mass distributions. 

The optimized selection criteria for the B Q — > 
J/^K s tt+tt- sample are p T {K%) > 1.0 GeV/c, ct(B) > 
30 pirn, and I > 0.50. This sample has less intrinsic back- 
ground than the B° — » J/tpK*°ir + ir~ sample since the 
invariant mass peak of the K s is narrower than that of 
the K*° and its decay vertex is additionally displaced 
from the J/tp decay vertex. To take advantage of this 
lower background, events that do not have enough SVX 
information to make a precise ct determination are in- 
cluded in a separate optimization which does not re- 
strict ct(B). The selection criteria for this sample are 
Pt(K^) > 1.9 GeV/c and I > 0.70. These two samples 
are combined in Fig. |2|(b). A fit yields 21.0 ± 6.3 signal 
candidates. 

The signal width is fixed in the fit to a = 11.3 MeV/c 2 , 
based upon the expected width from the Monte Carlo 
simulation scaled up by the ratio of widths between the 
data and Monte Carlo for the reference mode B° — » 
J/ipKg. Allowing this width to float results in a fitted 
width which is approximately half the expected width. 
The excess of events at the B° mass is robust across 
a wide range of selection criteria and is broader when 
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other selections are applied. The normalized mass dis- 
tribution (m B - 5.28 GeV/c 2 )/a mB has a = 0.67 ± 0.21 
and 16.2 ± 5.6 events. The unusual narrowness appears 
to be primarily an artifact of this particular set of selec- 
tion criteria which optimized the expected S 2 /(S+Bkg). 
Varying the width by ±20% affects the fitted signal size 
by less than 2%. 

The ratio of efficiencies, R e — 4.98, combined 
with 84.1 ± 9.9 B° — » J/ipK 1 ^ reference events 
leads to a ratio of branching ratios of BR(B° — > 
J/iPK s tt+tt~)/BR(B° -> J/i/>K%) = 1.24 ±0.40, where 
the error is statistical only. Using BR(B° -> J/ipK^) = 
(8.9 ±1.2) x 10" 4 leads to BR(B° -> J/tpK^n+ir ) = 
(1.1 ±0.4) x 10~ 3 . 

Unlike B° -> J/ipK*°7r+Tr- , B° -> J / ijj K s tt+ ir- 
shows evidence of several substructure contributions in 
addition to ip(2S)Kg candidates. The tt + tt~ and Kgir^ 1 
invariant mass plots shown in Fig. || have an excess 
of signal over background in the p° and K invari- 
ant mass regions, indicating possible contributions from 
B" -> J/ipK^p and B° -> J/ipK* + Tr~. The back- 
grounds are estimated from the tt + 7t _ and K^n^ invari- 
ant mass distributions of the candidates in the B mass 
sidebands. 
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FIG. 4: B° candidate invariant masses for events with a 
If candidate (a) and those with a p° candidate but no K* ± 
candidate (b). 



[4.1, 22.6] and [1.1, 15.6] for samples X and Y, respec- 
tively. Using the fitted number of signal events and the 
efficiencies for B° -> J/^K a s p a and B° -> J/^K*+ir- 
for each of the samples, the resulting branching ratios 
are BR(B° -> J/^K^p ) = (5.8 ± 3.1) x 10~ 4 and 



BR(B° -> J/^K*+tt 



.3 ±4.4) x 10~ 4 with a corre- 



lation coefficient of -0.43. These branching ratios assume 
that these two modes are the dominant contributions to 
the two samples and that they do not interfere in the 
overlap region of their K^tt + tt~ Dalitz plot. 

The J/ipir + ir~ invariant mass plot of Fig. ||(a) shows 4 
ip(2S) candidates on an expected background of 0.3. It is 
possible that B° — > J/tpK^p and J/i/jK* + tt~ come from 
B° -> J/ipKi(1270). Figure |(b) shows an excess in the 
if<j7r + 7r~~ invariant mass distribution near l£i(1270) but 
there is also a small excess of events at higher Kgir + ir~ 
invariant masses. The backgrounds are estimated from 
the J/^7r + 7r~ and Kgir + ir~ invariant mass distributions 
of the candidates in the B° invariant mass sidebands. 



FIG. 3: m(7r + 7r-) (a) and m^gvr*) (b) for J/ipK^iv + 'K~ 
events within ±2<r of the B° invariant mass. 

To fit for these contributions, the B° invariant mass 
peaks for two samples of events are considered. Events 
in sample X have a K*^ — > Kgir^ candidate with an in- 
variant mass within 0.892±0.051 GeV/c 2 . Sample Y con- 
tains events which have ap°-> 7r + 7r _ candidate with an 
invariant mass within 0.770±0.150 GeV/c 2 while exclud- 
ing events in sample X. The to(7t + 7t~) > 0.55 GeV/c 2 
requirement is not placed upon sample X. Neither sam- 
ple has any ip(2S) candidates. The invariant mass peaks 
of the B° candidates in these samples are fitted us- 
ing a Gaussian signal of fixed width and a linear back- 
ground. Figure ^ shows the results with sample X hav- 
ing 12.5 ± 4.6 fitted signal events and sample Y having 
8.5 ± 3.8 fitted signal events. 

Within ±2<7 of the B° mass, sample X has 21 events on 
a background of 9.2; sample Y has 14 events on a back- 
ground of 7.3. These numbers of events lead to Feldman- 
Cousins 95% confidence intervals Q for the signal size of 




3.63.8 4 4.24.44.64.8 5 

Mass [GeV/cl 



1 1.5 2 2.5 

Kg 3i* 7T Mass [GeV/c 2 ] 



FIG. 5: m(J/ipn + -K-) (a) and m(K^ir + ir~) (b) for 
J/ipKgiY + iT~ events within ±2cr of the B° invariant mass. 



The results are summarized in Table |. The domi- 
nant uncertainty in these branching ratios is the statisti- 
cal uncertainty due to the small signal size. Many of the 
systematic uncertainties cancel in the ratio of branching 
ratios with the reference mode. The systematic uncer- 
tainties that do not cleanly cancel are summarized in 
Table |n] and described as follows. 
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J/iPK* l >ir + TT- 


J/i>K u ir + ir- 


J/iPK u p u 


J/i>K* + TT- 


N obs 


85 


39 


14 


21 


Bkg 


54.0 


21.1 


7.32 


9.22 


c 

a fit 


oo.o ± y.y 


21.0 ±6.3 


8.5 ±3.8 


t n r 1 A r 
lz.0 ± 4.0 


Sfc 


[13.9, 50.8] 


[7.1, 31.8] 


[1.1, 15.6] 


[4.2, 22.6] 


$ref 


257 ± 18 


84.1 ±9.9 


84.1 ±9.9 


84.1 ±9.9 


Re 


3.75 


4.98 






BR 


(8.0 ±2.2 ±1.5) x 10" 4 


(1.1 ± 0.4 ± 0.2) x 10" 3 (5.8 ± 3.1 ± 1.2) x 10~ 4 


(8.3 ± 4.4 ± 1.7) x 10" 4 



TABLE I: Summary of results: Number of observed events (N b s ); Fitted background (Bkg); Fitted signal (£/«); Feldman- 
Cousins 95% confidence interval on the signal size (Sfc)', Fitted number of reference mode signal events (5 re /); Ratio of 
efficiencies (Re); and the branching ratio (BR) where the first uncertainty is statistical and the second is systematic. 



BR(B U -> J/4>Ktt + iv-) 


% Uncertainty 




Source of Uncertainty 


K = K*° 


K = K° 


Reference Mode BR 


11 


14 


Helicity Model 


9.9 


9.4 


Signal Width 


7.5 




Trigger Model 


5.0 


5.0 


Monte Carlo Composition 


5.0 


5.0 


B° Production Model 


2.5 


2.5 


Combined Uncertainties 


18 


18 


Without Ref. Mode BR 


15 


12 



TABLE II: Systematic uncertainties on BR(f?° — » 
J/^A-M'V+tt-)- 

The uncertainty on the reference mode branching ratio 
does not enter into the ratio of branching ratios but it 
is the dominant uncertainty for the branching ratio mea- 
surements. It is 11% for J/xPK*° and 14% for J/ipK° s . 
The signal modes have more helicity degrees of freedom 
than the reference modes and the relative contributions 
of possible helicity states are not known. This introduces 
an uncertainty in the efficiencies modeled with the Monte 
Carlo simulation of 9.9% for J/ipK*°Tr + w- and 9.4% for 
J/ipKgTr + Tr~ . There is a 5% uncertainty in both sig- 
nal modes due to uncertainties in the trigger model used 
in the Monte Carlo. Multiple decay modes could con- 
tribute to the final states studied here but they all have 
similar reconstruction efficiencies. Varying the relative 
compositions in the Monte Carlo results in a net un- 
certainty of 5%. The effect of differing pr(B) spectra 
from various B production models introduces a 3% uncer- 
tainty. Varying the fitted signal width and the width of 
the broad Gaussian shaped background by ±10% in the 
B° — > J/tpK*°Tr + ir~ sample results in a 7.5% variation in 
the fitted signal size and is thus included as a systematic 
uncertainty. Varying the signal width by ±20% in the 
B a — > J/tpKgir + ir~ sample has less than a 2% effect and 
thus is neglected. The J/xpKgp and J/xpK* + ir~ branch- 
ing ratios assume no interference in their overlap region 
in the Kgir + ir~ Dalitz plot. Completely constructive in- 
terference would increase their combined branching ratio 
by ~ 20%; a 10% systematic uncertainty is included in 



each mode to account for this possibility. 

The most statistically significant mode, J /ipK *°7r + 7r~, 
has a significance of 3.7a. J/ipKgTr + ir~ has a significance 
of 3.3ct; its submodes J/ipKgp and J/ipK* + n~ show 
hints of a signal but have less than 2a significance. The 
measured branching ratios are large enough that these 
modes should be visible in the data already recorded by 
the CLEO, Belle, and BaBar experiments. CDF should 
record hundreds of these decays in Run II. 
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